#include "simulation.h"

int main( int argc, char *argv[] )
{
    int iter;
    
    alea_init();
    
    for(iter=0;iter<16;++iter)
    {
        Simulation *sim = sim_create(2 + ( rand() % 500 ),
                                     2 + ( rand() % 500 ),
                                     5*alea(),
                                     6*alea());
        fprintf(stderr,
                "sim: %3dx%3d => [%g,%g]\n",
                sim->nx,
                sim->ny,
                sim->Lx,
                sim->Ly);
        
        field_set( sim->rho, 0.1 );
        sim_compute_pressure(sim);
        sim_compute_gradP(sim);
        sim_compute_momenta(sim);
        //sim_compute_velocities(sim);
        sim_compute_laplacian_rho(sim);
        sim_delete(sim);
    }
	return 0;
}
